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The observation of neutrinos from cosmic accelerators will be revolutionary. High energy neutrinos are closely 
connected to ultrahigh energy cosmic rays and their sources. Cosmic ray sources are likely to produce neutrinos 
and the propagation of ultrahigh cosmic rays from distant sources can generate PeV to ZeV neutrinos. We briefly 
review recent progress on the observations of ultrahigh energy cosmic rays and their implications for the future 
detections of high energy neutrinos. 



1. Introduction 

The detection of neutrinos from astrophysical 
sources will open a new window onto the work- 
ings of the highest energy cosmic accelerators and 
provide new ways to test neutrino interactions. 
Neutrino fluxes can be estimated if sources of ul- 
trahigh energy cosmic rays (UHECRs) are well 
understood. However, the sources of UHECRs, 
their location, cosmological evolution, and maxi- 
mum energy, as well as the injected composition, 
are presently unknown. 

Neutrinos can be generated at the sources of 
cosmic rays or during the propagation of UHE- 
CRs as they interact with the ambient pho- 
ton backgrounds jll2j . The neutrinos generated 
during UHECR propagation, named cosmogenic 
neutrinos, represent an almost guaranteed flux 
and have encouraged efforts to detect them for 
decades (see, e.g., [3 ). One important assump- 
tion for the existence of cosmogenic neutrinos, 
that cosmic rays are extragalactic at the high- 
est energies, has been verified by the detection 
of a feature consistent with the Greisen-Zatsepin- 
Kuzmin (GZK) cutoff '415' in the cosmic ray spec- 
trum |6|7) and by the indication of anisotropics 
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in the cosmic ray sky distribution at the high- 
est energies |8|9j . These findings herald a possi- 
ble resolution to the mystery behind the origin of 
UHECRs and the possibility of detecting ultra- 
high energy neutrinos in the near future. 

This optimistic view has been recently damped 
by the indication that UHECRs may be domi- 
nated by heavier nuclei [lOj. The cosmogenic 
neutrino fiux expected from heavy cosmic ray pri- 
maries can be much lower than if the primaries 
are protons at ultrahigh energies, making a de- 
tection extremely challenging for current obser- 
vatories. Conversely, if neutrinos are observed, 
they will test specific sets of cosmic ray source 
parameters. 

Below we review the current state of UHECR 
observations. In section 3, we describe the impact 
of these observations on predicted neutrino fluxes 
and conclude in section 4. 

2. Recent process in UHECRs 

The dominant component of cosmic rays ob- 
served on Earth are believed to originate in 
Galactic cosmic accelerators. A transition from 
Galactic to extragalactic cosmic rays should oc- 
cur somewhere between a PeV (= 10^^ eV) and 
a few EeV (= lO^^ eV). Above a few EeV, the 
so-called ultrahigh energy cosmic rays are most 
likely extragalactic. These are observed to reach 
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energies that exceed 10^*^ eV posing some inter- 
esting and challenging questions: Where do they 
come from? How can they be accelerated to 
such high energies? What kind of particles are 
they? What is the spatial distribution of their 
sources? What do they tell us about these ex- 
treme cosmic accelerators? How strong are the 
magnetic fields that they traverse on their way 
to Earth? How do they interact with the cosmic 
background radiation? What secondary particles 
are produced from these interactions? What can 
we learn about particle interactions at these oth- 
erwise inaccessible energies? Reviews on the sci- 
ence of ultrahigh energy cosmic rays can be found 
in [11112] , while for high energy neutrinos see [3] . 

Recent observations of UHECRs reveal a spec- 
trum whose shape supports the long-held no- 
tion that sources of UHECRs are extragalac- 
tic. As shown in Figure [T] the crucial spec- 
tral feature recently established at the highest 
energies is a steep decline in fiux above about 
30 EeV. This feature was first established by 
the HiRes observatory [6] and confirmed with 
higher statistics by the Pierre Auger Observa- 
tory [7] . The decline in flux is reminiscent of the 
effect of interactions between extragalactic cos- 
mic rays and the cosmic background radiation, 
named the Greisen-Zatsepin-Kuzmin (GZK) cut- 
off |4|5| . which causes cosmic ray protons above 
many tens of EeV to lose energy via pion pho- 
toproduction off cosmic backgrounds while cos- 
mic ray nuclei photodissociate. This feature was 
not seen in earlier observations with the AGASA 
array jlSI . Another important feature shown in 
Figure [T] is the hardening of the spectrum at a 
few EeV, called the ankle, which may be caused 
by the transition from Galactic to extragalactic 
cosmic rays or by propagation losses if UHECRs 
are mostly protons. 

Figure [T] also shows the observed spectrum fit 
by different models of UHECR sources (adapted 
from |llp. In the mixed composition and iron 
dominated models [15j , the ankle indicates a tran- 
sition from Galactic to extragalactic cosmic rays, 
the source evolution is similar to the star forma- 
tion rate (SFR), and the injection spectra are rel- 
atively hard (s '-^ 2 — 2.1). In the proton domi- 
nated models in the figure, the ankle is due to pair 
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Figure 1. Spectrum of UHECRs multiplied by 
observed by HiRes I |6 and Auger [M]. The dis- 
played error bars are statistical errors while the 
reported systematic error on the absolute energy 
scale is about 22%. Overlaid are simulated spec- 
tra obtained for different models of the Galactic 
to extragalactic transition and different injected 
chemical compositions and spectral indices, s. 



production propagation losses [T^, named "dip 
transition models" [T7], and the injection spec- 
tra are softer for a wide range of evolution mod- 
els. Models with proton primaries can also fit the 
spectrum with harder injection with a transition 
from Galactic to extragalactic at the ankle. 

The confirmed presence of a spectral feature 
similar to the predicted GZK cutoff, settles the 
question of whether acceleration in extragalac- 
tic sources can explain the high-energy spectrum, 
ending the need for exotic alternatives designed 
to avoid the GZK feature. However, the possibil- 
ity that the observed softening of the spectrum 
is mainly due to the maximum energy of accel- 
eration at the source, £'max, is not as easily dis- 
missed. A confirmation that the observed soft- 
ening is the GZK feature, awaits supporting evi- 
dence from the spectral shape, anisotropics, and 
composition and the observation of produced sec- 
ondaries such as neutrinos and photons. 
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The landmark measurement of a flux suppres- 
sion at the highest energies encourages the search 
for sources in the nearby extragalactic universe 
using the arrival directions of trans-GZK cos- 
mic rays (with energy above ~ 60 EeV). Above 
GZK energies, observable sources must lie within 
about 100 Mpc, the so-called GZK horizon or 
GZK sphere. At trans-GZK energies, light com- 
posite nuclei are promptly dissociated by cosmic 
background photons, while protons and iron nu- 
clei may reach us from sources at distances up 
to about 100 Mpc. Since matter is known to be 
distributed inhomogeneously within this distance 
scale, the cosmic ray arrival directions should ex- 
hibit an anisotropic distribution above the GZK 
energy threshold, provided intervening magnetic 
fields are not too strong. At the highest energies, 
the isotropic diffuse flux from sources far beyond 
this GZK horizon should be strongly suppressed. 

The most recent discussion of anisotropies in 
the sky distribution of UHECRs began with the 
report that the arrival directions of the 27 cosmic 
rays observed by Auger with energies above 57 
EeV exhibited a statistically significant correla- 
tion with the anisotropically distributed galaxies 
in the 12th VCV |I8j catalog of active galactic nu- 
clei ( AGN) |8l9j . The correlation was most signif- 
icant for AGN with redshifts z < 0.018 (distances 
< 75 Mpc) and within 3.1° separation angles. An 
independent dataset confirmed the anisotropy at 
a confidence level of over 99% |8I9| . The pre- 
scription established by the Auger collaboration 
tested the departure from isotropy given the VCV 
AGN coverage of the sky, not the hypothesis that 
the VCV AGN were the actual UHECR sources. 
A recent update of the anisotropy tests with 69 
events above 55 EeV (THj shows that the correla- 
tion with the VCV catalog is not as strong for the 
same parameters as the original period. The data 
after the prescription period shows a departure 
from isotropy at the Str level. With the currently 
estimated correlation fraction of 38%, a 5cr signif- 
icance will require at least another four years of 
Auger observations [19]. No corresponding corre- 
lation was observed in the northern hemisphere 
by HiRes [2U] where the distribution of their 13 
trans-GZK events is consistent with isotropy. 

The anisotropy reported by the test with the 



VCV catalog may indicate the effect of the large 
scale structure in the distribution of source har- 
boring galaxies or it may be due to a nearby 
source. An interesting possibility is the cluster 
of Auger events around the direction of Centau- 
rus A, the closest AGN (at - 3.8 Mpc). Only 
much higher statistics will tell if Cen A is the 
first UHECR source to be identified. 

The third key measurement that can help re- 
solve the mystery behind the origin of UHECRs 
is their composition as a function of energy ob- 
served on Earth. Composition measurements can 
be made directly up to energies of ^ 100 TeV 
with space-based experiments. For higher ener- 
gies, composition is derived from the observed 
development and particle content of the exten- 
sive airshower created by the primary cosmic ray 
when it interacts with the atmosphere. 

Observations of shower properties show the 
dominance of light nuclei around a few EeV. A 
surprising trend occurs in data by the Auger Ob- 
servatory above 10 EeV, a change toward heavy 
primaries is seen both in average position of the 
maximum of the showers as well as in the rms 
fiuctuations about the mean up to 40 EeV. As 
a mixture of different nuclei would increase the 
rms fluctuations, the observed narrow distribu- 
tion argues for a change toward a composition 
dominated by heavy nuclei. In contrast, the 
HiRes measurement of fluctuations remains closer 
to light primaries up to around 50 EeV. The two 
results are consistent within quoted errors, so the 
situation is currently unclear. 

Changes to hadronic interactions from current 
extrapolations provide a plausible alternative in- 
terpretation to the observed shower development 
behavior. Auger probes interactions above 100 
TeV center of mass, while hadronic interactions 
are only known around a TeV. The observation 
of anisotropies and secondary particles (neutrinos 
and gamma-rays) can lead to astrophysical con- 
straints on the composition of UHECRs, opening 
the possibility for the study of hadronic interac- 
tion cross sections, multiplicities, and other inter- 
action parameters at hundreds of TeV. 

The detailed composition of UHECRs is still 
to be understood, but it is clear that primaries 
are not dominated by photons |21l22j or neutrinos 
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|23)24j . Limits on the photon fraction place strin- 
gent hmits on models where UHECRs are gener- 
ated by the decay of super heavy dark matter and 
topological defects. Unfortunately, the uncertain- 
ties on the UHECR source composition, spec- 
trum, and redshift evolution translates to many 
orders of magnitude uncertainty in the expected 
cosmogenic neutrino flux as discussed next. 

3. Neutrinos from UHECRs 

Secondary neutrinos and photons can be pro- 
duced by UHECRs when they interact with ambi- 
ent baryonic matter and radiation fields inside the 
source or during their propagation from source to 
Earth. These particles travel in geodesies unaf- 
fected by magnetic fields and bear valuable in- 
formation of the birthplace of their progenitors. 
The quest for sources of UHECRs has thus long 
been associated with the detection of neutrinos 
and gamma rays that might pinpoint the position 
of the accelerators in the sky. 

The detection of these secondary particles is 
not straightforward however: first, the propaga- 
tion of gamma rays with energy exceeding several 
TeV is affected by their interaction with CMB 
and radio photons. These interactions lead to the 
production of high energy electron and positron 
pairs which in turn up-scatter CMB or radio 
photons by inverse Compton processes, initiating 
electromagnetic cascades. As a consequence, one 
does not expect to observe gamma rays of energy 
above ~ 100 TeV from sources located beyond a 
horizon of a few Mpc. Above EeV energies, pho- 
tons can again propagate over large distances, de- 
pending on the radio background, and can reach 
observable levels around tens of EeV. Secondary 
neutrinos are very useful because, unlike cosmic- 
rays and photons, they are not absorbed by the 
cosmic backgrounds while propagating through 
the Universe. In particular, they give a unique 
access to observing sources at PeV energies. How- 
ever, their small interaction cross-section makes it 
difficult to detect them on the Earth requiring the 
construction of km'^ or larger detectors. 

Neutrinos generated at UHECR sources have 
been investigated by a number of authors (be- 
ginning with [25126127] ). The normalization and 



the very existence of these secondaries highly de- 
pend on assumptions about the opacity of the ac- 
celeration region and on the shape of the injec- 
tion spectrum as well as on the phenomenolog- 
ical modeling of the acceleration. For instance, 
|27j obtain an estimate for the cosmic neutrino 
flux, by comparing the neutrino luminosity to the 
observed cosmic ray luminosity, in the specific 
case where the proton photo-meson optical depth 
equals unity. If the source is optically thick, [55] 
demonstrate that cosmic rays are not accelerated 
to the highest energies and neutrinos above E ~ 
EeV are sharply suppressed. 

The existence of secondary neutrinos from in- 
teractions during the propagation of cosmic rays 
is less uncertain, but it is also subject to large 
variations according to the injected spectral in- 
dex, chemical composition, maximum accelera- 
tion energy, and source evolution history. A 
number of authors have estimated the cosmo- 
genic neutrino flux with varying assumptions 
(e.g., |29iaOI81l32l83l84i:i5i:ffil87] 'l. Figure[2] sum- 
marizes the effects of different assumptions about 
the UHECR source evolution, the Galactic to ex- 
tragalactic transition, the injected chemical com- 
position, and i?max7 on the cosmogenic neutrino 
flux. It demonstrates that the parameter space is 
currently poorly constrained with uncertainties of 
several orders of magnitude in the predicted flux. 
UHECR models with large proton -Eniax(> 100 
EeV), source evolution corresponding to the star 
formation history or the GRB rate evolution, dip 
or ankle transition models, and pure proton or 
mixed 'Galactic' compositions are shaded in grey 
in Figure [2] and give detectable fluxes in the EeV 
range with 0.06—0.2 neutrino per year at IceCube 
and 0.03 — 0.06 neutrino per year for the Auger 
Observatory. If EeV neutrinos are detected, PeV 
information can help select between competing 
models of cosmic ray composition at the highest 
energy and the Galactic to extragalactic transi- 
tion at ankle energies. With improved sensitivity, 
ZeV (=10^^ eV) neutrino observatories, such as 
ANITA and JEM-EUSO could explore the maxi- 
mum acceleration energy. 

Due to the delay induced by EGMF on charged 
cosmic rays, secondary neutrinos and photons 
should not be detected in time coincidence with 
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Figure 2. Cosmogenic neutrino flux for all fla- 
vors, for different UHECR parameters compared 
to instrument sensitivities. Dash-dot line corre- 
sponds to a strong source evolution case with a 
pure proton composition, dip transition model, 
and -Emax — 3 ZeV. Uniform source evolution 
with iron rich (30%) composition and -Ez.max < ^ 
10 EeV is shown in the dotted line and dashed line 
represents pure iron injection and -Bz,max = Z 
100 EeV. Grey shaded range brackets dip and 
ankle transition models, with evolution of star 
formation history for z < 4, pure proton and 
mixed 'Galactic' compositions, and large proton 
E^max(> 100 EeV)). Including the uniform source 
evolution would broaden the shaded area down 
to the solid line below it. Current experimental 
limits (labeled lines on top) assume 90% confi- 
dence level and full mixing neutrino oscillation. 
The differential limit and the integral flux limit 
on a pure spectrum (straight lines) are pre- 
sented for IceCube 22 lines [38], ANITA-II [39] 
and Auger South |23j . For future instruments, 
we present the projected instrument sensitivities 
(dashed lines) for IceCube 80 lines (acceptances 
from S. Yoshida, private communication, see also 
[40]), and for JEM-EUSO j41j. 



UHECRs if the sources are not continuously emit- 
ting particles, but are transient such as GRBs and 
young magnetars. 

4. Conclusion 

The possibility of observing high energy neutri- 
nos is intimately related to the resolution of the 
long standing mystery of the origin of UHECRs. 
To discover the origin of UHECRs will require 
a coordinated approach on three complementary 
fronts: the direct UHECR frontier, the transition 
region between the knee and the ankle, and the 
multi-messenger interface with high-energy pho- 
tons and neutrinos. 

Current data suggest that watershed 
anisotropics will only become clear above 60 
EeV and that very large statistics with good 
angular and energy resolution will be required. 
The Auger Observatory (located in Mendoza, 
Argentina), will add 7 x Vd^krv? .sr each year of 
exposure to the southern sky, while the Telescope 
Array (located in Utah, USA) will add about 
2 X V^^krr? .sr each year in the North. Current 
technologies can reach a goal of another order 
of magnitude if deployed by bold scientists over 
very large areas (e.g.. Auger North). New tech- 
nologies may ease the need for large number of 
detector units to cover similarly large areas. Fu- 
ture space observatories (e.g., JEM-EUSO, OWL, 
Super-EUSO) promise a new avenue to reach the 
necessary high statistics especially if improved 
photon detection technologies are achieved. 

Existing and upcoming high energy neutrino 
detectors roughly cover three energy ranges: PeV 
(= 1015 eV), EeV (= lO^^ eV), and ZeV (= 
10^1 eV). ANTARES, IceCube, and the future 
KM3Net are large water or ice cubic detectors 
that aim at observing events around PeV ener- 
gies. IceCube will also have a very good sensitiv- 
ity at higher energies, and will ultimately be able 
to cover a wide energy range from about 1 PeV 
to ~ 10 EeV. Experiments primarily dedicated 
to the detection of cosmic rays like the Pierre 
Auger Observatory and the Telescope Array have 
their best neutrino sensitivities in the EeV energy 
range. The radio telescope ANITA and the fluo- 
rescence telescope JEM-EUSO are most effective 
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at the highest energy neutrinos around 0.1 ZeV. 

With a significant increase in the integrated ex- 
posure to neutrinos and cosmic rays, next gener- 
ation observatories may reach the sensitivity nec- 
essary to achieve charged particle astronomy and 
to observe ultrahigh energy photons and neutri- 
nos, which will further illuminate the workings of 
the universe at the most extreme energies. 

REFERENCES 

1. V. S. Berezinsky, G. T. Zatsepin, Physics Let- 
ters B 28 (1969) 423. 

2. F. W. Stecker, ApJ 228 (1979) 919. 

3. L. A. Anchordoqui, T. Montaruli, Ann. Rev. 
Nucl. Part. Sci. 60, 129. 

4. K. Greisen, Phys. Rev. Lett. 16 (1966) 748. 

5. G. Zatsepin, V. Kuzmin, J. Exp. Thcor. Phys. 
Lett. 4 (1966) 78. 

6. R. U. Abbasi, et al., Astroparticle Physics 32 

(2009) 53. 

7. J. Abraham, et al., Phys. Rev. Lett. 101 
(2008) 061101. 

8. J. Abraham, et al.. Science 318 (2007) 938. 

9. J. Abraham, et al., Astroparticle Physics 29 
(2008) 188. 

10. J. Abraham, et al., Phys. Rev. Lett. 104 

(2010) 091101. 

11. K. Kotera, A. V. Olinto, Annu. Rev. Astron. 
Astrophys. v. 49 (2011), |arXiv:1101.42 56vl 

12. A. Letessier-Selvon, T. Stanev, Rev. Mod. 
Phys. (2011). 

13. M. Takeda, et al.. Physical Review Letters 81 
(1998) 1163. 

14. J. Abraham, et al., Physics Letters B 685 
(2010) 239. 

15. D. AUard, E. Parizot, A. V. Ohnto, Astropar- 
ticle Physics 27 (2007) 61. 

16. V. S. Berezinsky, S. I. Grigorieva, A&A 199 
(1988) 1. 

17. V. Berezinsky, A. Gazizov, S. Grigorieva, 
Phys. Rev. D 74 (4) (2006) 043005 

18. M.-P. Veron-Cetty, P. Veron, A&A 455 
(2006) 773. 

19. P. Abreu, et al., Astroparticle Physics 34 
(2010) 314. 

20. R. U. Abbasi, et al., Astroparticle Physics 30 
(2008) 175. 



21. J. Abraham, et al., Astropart. Phys. 29 (2008) 
243. 

22. J. Abraham, et al, Astropart. Phys. 31 (2009) 
399. 

23. J. Abraham, et al., Phys. Rev. D 79 (10) 
(2009) 102001. 

24. R. U. Abbasi, et al., Astrophys. J. 684 (2008) 
790. 

25. A. P. Szabo, R. J. Protheroe, Astroparticle 
Physics 2 (1994) 375. 

26. J. P. Rachen, P. Meszaros, Phys. Rev. D 
58 (12) (1998) 123005 

27. E. Waxman, J. BahcaU, Phys. Rev. D 59 (2) 
(1999) 023002 

28. D. AUard, R. J. Protheroe, A&A 502 (3) 

(2009) 803. 

29. R. Engel, D. Seckel, T. Stanev, Phys. Rev. D 
64 (9) (2001) 093010 

30. M. Ave, et al., Astropart. Phys. 23 (2005) 19. 

31. D. Seckel, T. Stanev, Physical Review Letters 
95 (14) (2005) 141101 

32. D. Hooper, A. Taylor, S. Sarkar, Astroparticle 
Physics 23 (2005) 11. 

33. V. Berezinsky, Nuclear Physics B Proceedings 
Supplements 151 (2006) 260. 

34. T. Stanev, D. de Marco, M. A. Malkan, F. W. 
Stecker, Phys. Rev. D 73 (4) (2006) 043003. 

35. D. AUard, M. Ave, N. Busca, M. A. Malkan, 
A. V. Olinto, E. Parizot, F. W. Stecker, 
T. Yamamoto, JCAP 9 (2006) 5. 

36. H. Takami, K. Murase, S. Nagataki, K. Sato, 
Astroparticle Physics 31 (2009) 201. 

37. K. Kotera, D. AUard, A. V. Ohnto, JCAP 10 

(2010) 13. 

38. R. U. Abbasi, et al.. Physical Review Letters 
104 (16) (2010) 16110L 

39. P. W. Gorham, et al., 'arXiv: 1003.2961] 

40. A. Karle, IceCube, arXiv: 1003.5715 

41. Y. Takahashi, et al.. The Jem-Euso Mission, 
New J. Phys. 11 (2009) 065009. 



